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Research on graphene and other two-dimensional (2D) materials, such as silicene, germanene, 
phosphorene, hexagonal boron nitride (h-BN), graphitic carbon nitride (g-C 3 N 4 ), graphitic zinc ox¬ 
ide (g-ZnO) and molybdenum disulphide (M 0 S 2 ), has recently received considerable interest owing 
to their outstanding optoelectronic properties and wide applications. Looking beyond this field, 
combining the electronic structures of 2D materials in ultrathin van der Waals heterojunctions 
has also emerged to widely study theoretically and experimentally to explore some new proper¬ 
ties and potential applications beyond their single components. Here, this article reviews our re¬ 
cent theoretical studies on the structural, electronic, electrical and optical properties of 2D van der 
Waals heterojunctions using density functional theory calculations, including the Graphene/Silicene, 
Graphene/Phosphorene, Graphene/g-ZnO, Graphene/MoS 2 and g-C3N4/MoS2 heterojunctions. 
Our theoretical simulations, designs and calculations show that novel 2D van der Waals hetero¬ 
junctions provide a promising future for electronic, electrochemical, photovoltaic, photoresponsive 
and memory devices in the experiments. 


INTRODUCTION 

Two-dimensional (2D) ultrathin materials, such as 
graphene [1-3], silicene [4-7], germanene [8-10], phos¬ 
phorene [11-13], hexagonal boron nitride (h-BN) [14- 
16], graphitic carbon nitride (g-C 3 N 4 ) [17-19], graphitic 
zinc oxide (g-ZnO) [20-22] and molybdenum disulphide 
(M 0 S 2 ) [23-25], have received considerable interest re¬ 
cently owing to their outstanding properties and wide 
applications. There already have been many review ar¬ 
ticles [26-32] for the research on 2D materials over the 
past several years. 

Graphene [2], a 2D sp^-hybridized carbon monolayer, 
is known to have remarkable electronic properties, such 
as a high carrier mobility, but the absence of a bandgap 
limits its applications of large-off current and high on-off 
ratio for graphene-based electronic devices. The same 
limitation also exists in silicene [4] and germanene [8], 
which have most similar remarkable electronic proper¬ 
ties to graphene but with buckled honeycomb structures. 
On the other hand, most of graphene’s derivatives, such 
as graphene oxide [33-35], graphane [36-38] and fluoro- 
graphene [39-41] are all semiconductors but with lower 
carrier mobility, and have also been studied experimen¬ 
tally and theoretically. 

Phosphorene, a new two dimensional (2D) elemental 
monolayer [11], has recently been experimentally isolated 
through mechanical exfoliation from bulk black phospho¬ 
rus. Phosphorene exhibits some remarkable electronic 
properties superior to graphene. For example, phospho¬ 
rene is a direct semiconductor with a high hole mobil¬ 
ity [12], showing the drain current modulation up to 10® 
in nanoelectronics [13]. Besides graphene, phosphorene is 


the only stable elemental 2D monolayer which can be me¬ 
chanically exfoliated experimentally [12]. Most recently, 
arsenene and antimonene have also predicted theoreti¬ 
cally [42]. Monolayer arsenene and antimonene are in¬ 
direct wide-band-gap semiconductors, but become direct 
band-gap semiconductors under strain. Owing to these 
band-gap transitions, they can be used in nanoelectronic 
and optoelectronic devices. 

2D h-BN monolayer [14], so-called ’’white graphene”, 
is comprised of alternating boron and nitrogen atoms in 
a sp^-bonded 2D honeycomb arrangement and isolated 
from bulk BN. It has a wide bandgap used for its ec- 
cellent electrical insulations, high thermal conductivity 
and superior lubricant properties, and can be useful as 
an important complementary 2D dielectric substrate for 
graphene electronics. 

2D g-C 3 N 4 monolayer [17], which can be made by poly¬ 
merization of cyanamide, dicyandiamide or melamine, is 
a semiconductor with a bandgap of about 2.7 eV, as 
a new polymeric photocatalyst for hydrogen production 
from water splitting under visible light, which has been 
demonstrated to be a be a potential candidate for solar 
cell absorber and photovoltaic materials. 

ZnO nanofilms and nanoparticles have been widely 
used with graphene as composite materials for light- 
emitting diodes [43], photocatalysts [44], quantum 
dots [45] and field effect transistors [46]. However, first- 
principles calculations [20, 21] have predicted and then 
confirmed experimentally [22] that the ZnO(OOOl) film 
prefers a graphitic honeycomb structure as the layer num¬ 
ber is reduced. Furthermore, 2D g-ZnO monolayer itself 
is also semiconducting with a wide band gap of 3.57 eV, 
showing many interesting electronic and magnetic prop¬ 
erties [47-49]. 
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2D M 0 S 2 is an important kind of transition metal 
dichalcogenides (TMDCs), which has been widely stud¬ 
ied experimentally and theoretically [23-25]. M 0 S 2 has a 
direct band gap of about 2 eV, showing a high on-off cur¬ 
rent ratio with a high carrier mobility of around 200-500 
cm^/(Vs) in nanotransistors. Most of other TMDCs [29], 
such as MoSe 2 , WS 2 and WSe 2 , but have attracted rel¬ 
atively little attention despite the fact that they have 
similar electronic structures to M 0 S 2 . 

In parallel with the great efforts on 2D materials, an¬ 
other research field has recently emerged has been gain¬ 
ing strength over the past several years. It deals with 
ultrathin 2D van der Waals heterojunctions and devices 
made layer by layer from 2D materials, stacking differ¬ 
ent 2D monolayers on the top of each other. Strong co¬ 
valent bonds provide in-plane stability of 2D monolay¬ 
ers, whereas relatively weak van der Waals forces are 
sufficient to keep the stack together. The possibility 
of making 2D van der Waals heterojunctions has been 
demonstrated experimentally for electronic, electrochem¬ 
ical, photovoltaic, photoresponsive and memory devices. 
Recently, there only have been several review articles [50- 
53] involved for the research on 2D van der Waals het- 
ero junctions. 

Up to now, significant experimental and theoreti¬ 
cal efforts have been made to characterize and ex¬ 
plore new structures, properties and applications of 
2D van der Waals heterojunctions, especially graphene 
based heterojunctions, such as graphene/hexagonal 
boron nitride (G/h-BN) [54-58], graphene/silicene 
(G/S) [59-61], graphene/phosphorene (G/P) [62-64], 
graphene/graphitic carbon nitride (G/g-G 3 N 4 ) [65- 
67], graphene/graphitic zinc oxide (G/g-ZnO) [68-70], 
graphene/molybdenum disulphide (G/M 0 S 2 ) [71-75], 
graphene/molybdenum disulfide (G/MoSe 2 ) [76-78], sil- 
icene/hexagonal boron nitride (S/h-BN) [79-81], sil- 
icene/molybdenum disulphide (S/M 0 S 2 ) [82-84], and 
TMDGs based hetero junctions [85-87]. These van der 
Waals heterojunctions show much more new properties 
far beyond their single components. Furthermore, most 
of them are ideal substrates for each other to preserve 
their intrinsic electronic properties due to weak van der 
Waals interactions at the interface. 

Hybrid G/h-BN heterojunctions [54-58] are the most 
widely studied topic in 2D van der Waals heterojunctions 
during the past few years. Theoretically, Giovannetti et 
al. [54] proved that substrate-induced band gaps can be 
opened at the Dirac point of graphene on h-BN surfaces 
and the gap values are tunable as varying the the in¬ 
terfacial distance, which greatly improves the room tem¬ 
perature pinch-off characteristics of graphene-based field 
effect transistors. Experimentally, Dean et al. [55] first in 
2010 introduced atomically thin 2D van der Waals het¬ 
erojunctions in which graphene is on top of h-BN layers, 
and shown that h-BN can be used as an excellent comple¬ 
mentary 2D dielectric substrate for graphene electronics. 


In the present work, we review our recent theoretical 
studies on the structural, electronic, optical and electrical 
properties of 2D van der Waals hetero junctions with den¬ 
sity functional theory calculations and shows that these 
2D van der Waals heterojunctions provide a promising 
future for electronic, electrochemical, photovoltaic, pho¬ 
toresponsive and memory devices [50-53]. 

THEORETICAL MODELS 

In order to simulate hybrid composite structures of 2D 
van der Waals heterojunctions made layer by layer, we 
need to choose two special supercells for each 2D mate¬ 
rial with different crystal systems and lattice constants, 
making them have a smaller lattice mismatch for each 
other. We have written a small program named Lat- 
ticeMatch(See the supplementary material) to construct 
theoretical models for 2D van der Waals heterojunctions 
with a smaller lattice mismatch for different 2D mate¬ 
rials. Its basic idea is based on the rotation matrix of 
extended orthogonal lattices for different 2D materials 
to search their minimum matching supercell models in 
2D van der Waals heterojunctions. 

COMPUTATIONAL METHODS 

Our first-principles calculations are based on the den¬ 
sity functional theory (DFT) implemented in the VASP 
package. [88] All the geometry structures are fully relaxed 
until energy and forces are converged to I0“® eV and O.OI 
eV /A, respectively. Dipole correction is employed to can¬ 
cel the errors of electrostatic potential, atomic forces and 
total energy, caused by periodic boundary condition [89]. 
Very fine points in the surface Brillouin zone sampled 
with a regular mesh are used for calculating the tiny band 
gaps at the Dirac points of graphene and silicene in our 
recently studied 2D van der Waals heterojunctions. 

Because conventional exchange-correlation functionals 
in the DFT calculations, such as the generalized gra¬ 
dient approximation of Perdew, Burke, and Ernzerhof 
(GGA-PBE) [90], fail to adequately describe the long- 
range weak van der Waals interactions, different meth¬ 
ods, such as semi-empirical long-range dispersion correc¬ 
tion (DFT-D) proposed by Grimme et al. [91-93] and 
non-empirical van der Waals density functional (vdW- 
DE) scheme proposed by Dion et al. [94-96], have been 
developed. Table I shows different methods implemented 
in VASP calculated equilibrium interfacial distance and 
binding energy per carbon atom for bilayer graphene (a 
= b = 2.47 A) compared with experimental measure¬ 
ments [97, 98]. In our recent theoretical studies on 2D van 
der Waals heterojunctions, we choose the semi-empirical 
DPT-D2 method [92] due to its widely used and good de¬ 
scription of long-range weak van der Waals interactions in 
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TABLE I: Different methods implemented in VASP calculated 
equilibrium interfacial distance Dq (A) and binding energy 
per carbon atom Et (meV) for bilayer graphene compared 
with experimental measurements. 


Bilayer graphene 

Do 

Et 

GGA (PBE) 

3.90 

-0.04 

DFT-D2 (PBE) 

3.25 

-25.2 

vdW-DF (PBE) 

3.47 

-33.6 

vdW-DF (optB88) 

3.32 

-32.0 

Exp. [97, 98] 

3.35 

-26.0 


TABLE II: DFT-D2 calculated equilibrium interfacial dis¬ 
tance Do (A) and binding energy per carbon atom Ei, (meV) 
in our studied 2D van der Waals heterojunctions 


DFT-D2 

Do 

Ei 

Graphene / Silicene 

3.37 

-24.7 

Graphene/Phosphorene 

3.43 

-24.7 

Graphene/g-ZnO 

3.14 

-51.0 

Graphene/MoS 2 

3.37 

-20.5 

g-G3N4/MoS2 

2.97 

-22.4 


molecular surface adsorption and layered structure sys¬ 
tems [99-109]. 

The charge carrier (hole and electron) concentration of 
doping graphene in 2D van der Waals heterojunctions can 
be estimated by the linear dispersion around the Dirac 
point of graphene [3] 


-^h/e = 


TT{hVFY 


( 1 ) 


where AEf, is the shift of graphene’s Dirac point {Eo) 
relative to the Fermi level (Ep), that is AE]j = Ejj - 
Ep. 

To study the optical properties of 2D van der Waals 
heteroj unctions, the frequency-dependent dielectric ma¬ 
trix is calculated [110]. The imaginary part of dielectric 
matrix is determined by a summation over states as 


ea/3 4 II 2u;k<5(eek - e.k - w) 

C,D,k 

^ (Mck+e^ q l/^uk) (Mck+e^ q\f^vU.) 


( 2 ) 


where, D is the volume of the primitive cell, Wk is the 
k point weight, c and v are the conduction and valence 
band states respectively, Cck and fick are the eigenvalues 
and wavefunctions at the k point respectively, and are 
the unit vectors for the three Cartesian directions. In 
order to accurately calculate the optical properties of 2D 
van der Waals heterojunctions, a large four times regular 
mesh for the surface Brillouin zone, a large number of 
empty conduction band states (two times more than the 
number of valence band) and frequency grid points are 
adopted in optical properties calculations compared with 
conventional electronic structure calculations. 


VAN DER WAALS HETEROJUNCTIONS 

We hrst check some common properties in our 
recently studied 2D van der Waals heterojunc¬ 
tions, including the Graphene/Silicene (G/S) [59], 
Graphene/Phosphorene (G/P) [62], Graphene/g-ZnO 


(G/g-ZnO) [68], Graphene/MoS 2 (G/M 0 S 2 ) [75] 

and g-G 3 N 4 /MoS 2 [87] heterojunctions as shown in 
Fig. 1, 2, 3, 4 and 5, respectively. We find that the 
weak van der Waals interlayer interactions in 2D van der 
Waals heterojunctions can induce new properties and 
phenomena, such as bandgap opening, charge transfer 
and new optical absorption. 

Firstly, we Hnd that weak van der Waals interactions 
domain at the interfaces of these 2D van der Waals het¬ 
erojunctions. Table II shows our recently studied 2D van 
der Waals heteroj unctions all have typical vdW equilib¬ 
rium spacings with corresponding small binding energies, 
thus, they can be used as ideal substrates for each other 
with their intrinsic electronic structures undisturbed in 
2D van der Waals heterojunctions. Furthermore, we con¬ 
firm that the small lattice mismatch of about 2% for spe¬ 
cial supercells of 2D materials obtained by our Lattice- 
Match program has little effect on their electronic proper¬ 
ties in these 2D van der Waals heterojunctions as shown 
their electronic band structures and density of states in 
Fig. lb, 2b, 3b, 4b and 5b. 

Secondly, small band gaps are opened at the Dirac 
points of graphene in graphene based 2D van der Waals 
heteroj unctions due to weak van der Waals interactions 
at the interfaces, though its linear Dirac-like dispersion 
relation around the Fermi level is still preserved. Notice 
that induced graphene band gaps in 2D van der Waals 
heteroj unctions are typically sensitive to other external 
conditions, such as interlayer separation [56], showing 
that the band gap values increase gradually with the in¬ 
terlayer separation decrease, thus tunable, with potential 
applications for graphene-based Held effect transistors. 

Thirdly, charge redistribution and transfer always oc¬ 
cur at the interfaces in 2D van der Waals heterojunc¬ 
tions due to the interlayer coupling. In our recently 
studied 2D van der Waals heteroj unctions, the G/S het¬ 
erojunction belongs to the Metal/Metal type, the G/g- 
ZnO, G/P, and G/M 0 S 2 heterojunctions belong to the 
Metal/Semiconductor type, the g-G 3 N 4 /MoS 2 hetero¬ 
junction belongs to the Semiconductor/Semiconductor 
type. For metal, they have different work functions (The 
energy difference between the vacuum level and the Fermi 
level). For semiconductors, they have different ioniza- 
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FIG. 1: (Color online) (a) Geometric structures of G/S heterojunction (top and side views). The gray and yellow balls denote 
carbon and silicon atoms, respectively, (b) Electronic band structures of silicene, graphene and G/S heterojunction. The Fermi 
level is set to zero and marked by green dotted lines, (c) The doping charge carrier concentrations A^h/e (10^^ cm~^) of silicene 
and graphene of in G/S heterojunction as a function of interfacial distance D (A). The binding energies Eb (meV) per atom 
and the Dirac points AEo (eV) of silicene and graphene shift relative to the Fermi level are shown in the inset. The equilibrium 
spacing is denoted by green dotted lines, (d) XY-averaged electrostatic potentials of G/S heterojunction at different interfacial 
distances D (A) in the Z direction. Depths of potential wells of graphene and silicene are shown in the inset, (e) Imaginary 
part (e") of frequency {E = hui) dependent dielectric function (parallel and perpendicular) for pristine silicene and graphene 
monolayers as well as corresponding G/S heterojunction. (Reproduced from Ref. [59], with permission from American Institute 
of Physics.) 


tions (The energy difference between the vacuum level 
and the conduction band minimum) and nucleophilics 
(The energy difference between the vacuum level and 
the valence band maximum) potentials. Based on the 
Schottky-Mott model [111], different work functions, ion¬ 
izations and nucleophilics can induce charge transfer as 
well as well-separated electron-hole pairs at the inter¬ 
faces in 2D van der Waals heterojunctions (Fig. Ic, 2e, 
3c and 5a). Furtheremore, inhomogeneous surfaces also 
can results in charge redistribution and form intralayer 
electron-hole puddles at the interfaces in 2D van der 
Waals heterojunctions as shown in Fig. 3b. These new 
properties of charge redistribution and transfer at the in¬ 
terfaces can significantly enhance the conductivity and 
generate new catalytic activities in 2D van der Waals 
heterojunctions. 

Fourthly, interlayer interaction [19] and charge trans¬ 
fer [67] in 2D van der Waals heterojunctions may in¬ 
duce new optical transitions, though pristine monolayers 
themselves display outstanding optical properties. We 
find that 2D van der Waals hetero junctions always ex¬ 


hibit wider absorption range and stronger optical absorp¬ 
tion compared with their single monolayers as shown in 
Fig. le, 2e, 3d, 4d and 5c, because electrons can now be 
directly excited between different layers at the interfaces 
of heterojunctions. These new optical transitions in 2D 
van der Waals heterojunctions are expected to be with 
a great potential in photocatalytic and photovoltaic de¬ 
vices. 

In the following of this review, we primarily focus on 
some special features of our recently studied 2D van der 
Waals heteroj unctions that distinguish them from other 
previously studied 2D van der Waals heterojunctions. 


Graphene / Silicene 

Graphene and silicene are metallic but they have dif¬ 
ferent work functions. Therefore, both silicene and 
graphene can be doped in hybrid G/S heterojunction 
as plotted in Fig. lb, and their Dirac points respec¬ 
tively shift 0.12 eV above and below the Fermi lev- 
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FIG. 2: (Color online) (a) Geometric structures of G/P heterojunction (top and side views). The gray and violet balls denote 
carbon and phosphorus atoms, respectively, (b) Electronic band structures of graphene, phosphorene and G/P heteroj unction. 
The Fermi level is set to zero and marked by green dotted lines, (c) Electronic band structures of G/P heteroj unction at different 
interfacial distances D = 2.8, 3.0, 3.5, 4.0 and 4.5 A. (d) Schottky barriers ^Bn, ^Bp and ^Bn+^Bp in G/P heteroj unction as 
a function of interfacial distance, (e) Imaginary part of frequency dependent dielectric function (parallel and perpendicular) 
for pristine graphene and phosphorene monolayers as well as corresponding G/P heterojunction. Differential charge density 
(0.002 e/K?) of G/P heterojunction is shown in the insert. The red and blue regions indicate electron increase and decrease, 
respectively. (Reproduced from Ref. [62], with permission from Royal Society of Chemistry.) 


els in hybrid G/S heterojunction, inducing weak p-type 
and n-type doping of silicene and graphene, respectively. 
Our calculated charge carrier concentrations are fV/j(S) 
= 3.8x10^^ and A/.(G) = 1.6x10^^ for sil¬ 

icene and graphene in hybrid S/G heterojunction, respec¬ 
tively. Notice that these charge carrier concentrations 
are more than 2 orders of magnitude larger than the in¬ 
trinsic charge carrier concentration of graphene at room 
temperature (n = irk^T^/6hiyp = 6x10^° cm~^). Fur¬ 
thermore, we find that the charge carrier concentrations 
of both silicene and graphene in hybrid G/S heterojunc¬ 
tion can be tuned via the interfacial spacing as shown 
in Fig. Ic. Their charge carrier concentrations decrease 
gradually with the interlayer separation increase and a 
conversion of doping types of silicene and graphene oc¬ 
curs when their interfacial distance artificially increases 
to above 4.2 A. Therefore, the self-doping phenomenon 
in hybrid G/S heteroj unction provides a effective and 
tunable way for new optoelectronic devices. 

Here, we reveal the origin intrinsic mechanism of 


tunable self-doping in hybrid G/S heteroj unction. We 
observe that AEo is already unchanged for both sil¬ 
icene and graphene as the interfacial distance increases 
larger than 4.6 A, but AEb converges to different values 
{AEd{S) = -0.01 eF and AEd{G) = 0.01 eV). It can 
be interpreted by their different work functions, Wf{S) 
= 4.6 eV and Wf{G) = 4.4 eV. Based on the Schottky- 
Mott model [111], few electrons should be transferred 
from graphene to silicene, resulting in weak n-type and p- 
type doping of silicene and graphene, respectively. When 
the interfacial distance decreases, weak overlap of elec¬ 
tronic states between silicene and graphene is enhanced, 
which increase the charge transfer and electric double 
layer at the G/S interface, shifting up and down the en¬ 
ergy levels of silicene and graphene respectively as shown 
in Fig. Id. Furthermore, the tunneling energy barrier 
for electrons at the G/S interface is also reduced by 
their interface interactions. Therefore, when silicene and 
graphene are close to each other (smaller than 4.2 A), 
few electrons are transferred from silicene to graphene. 
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FIG. 3: (Color online) (a) Geometric structures structures of G/g-ZnO heterojunction (top and side views). The gray, red 
and blue balls denote carbon, oxygen and zinc atoms, respectively, (b) Electronic band structures of G/g-ZnO, G/Al doped 
g-ZnO, and G/Li doped g-ZnO heterojunctions. The Fermi level is set to zero and marked by red dotted lines, (c) Differential 
charge density (0.002 e/A®) for G/g-ZnO, G/Al doped g-ZnO, and G/Li doped g-ZnO heteroj unctions. The red and blue 
regions indicate electron increase and decrease, respectively, (d) Imaginary part (e") of dielectric function (perpendicular) 
for graphene, g-ZnO monolayers and G/g-ZnO heteroj unctions. (Reproduced from Ref. [68], with permission from American 
Institute of Physics.) 


resulting in weak p-type and n-type doping of silicene 
and graphene, respectively. That is why a reversion of 
doping types of silicene and graphene occurs as the in¬ 
terfacial distance is 4.2 A. In previous studies, similar 
feature of tunable doping types and charge carrier con¬ 
centration of graphene on metal substrates [112] have 
been also reported experimentally and theoretically. 


Graphene/Phosphorene 

Monolayer phosphorene is semiconducting with a di¬ 
rect band gap of 0.85 eV (Fig. 2b), which is different form 
graphene. Therefore, a Schottky contact can be formed 
between metallic graphene and semiconducting phospho¬ 
rene as shown in Fig. 2d. Based on the Schottky-Mott 
model [111] at the metal/semiconductor interface [113], 
a n-type Schottky barrier ($b„) is defined as the energy 
difference between the Fermi level (Ep) and the conduc- 
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FIG. 4: (Color online) (a) Geometric strnctures of G/M 0 S 2 and G/M 0 S 2 /G heterojunctions. The gray, yellow and blue balls 
denote carbon, sulfur and molybdenum atoms, respectively, (b) Electronic band structures of graphene and M 0 S 2 as well as 
corresponding G/M 0 S 2 and G/M 0 S 2 /G heterojunctions. The Fermi level is set to zero and marked by green dotted lines, 
(c) The doping charge carrier concentrations (10^^ cm~^) of graphene in G/M 0 S 2 and G/M 0 S 2 /G heterojunctions as a 
function of vertical electric field E (I//A). The energy shift AEd (eV) of graphene’s Dirac point relative to the Fermi level is 
shown in the inset, (d) Imaginary part of dielectric function {e") of pristine graphene M 0 S 2 monolayers as well as corresponding 
G/M 0 S 2 and G/M 0 S 2 /G heteroj unctions. (Reproduced from Ref. [75]) 


tion band minimum {Eq), that is ^Bn = Eq - Ep. Sim¬ 
ilarly, a p-type Schottky barrier {^Bp) is defined as the 
energy difference between the Fermi level {Ep) and va¬ 
lence band maximum (Ey), that is = Ep - Ey. 
Furthermore, the sum of two types (n- and p-type) of 
Schottky barrier is approximately equal to the band gap 
value {Eg) of semiconductor, that is ^Bn + ^Bp ^ Eg- 
As shown in Fig. 2c, the Dirac point of graphene moves 
from the conduction band to the valence band of phos- 
phorene as the interfacial distance decreases from 4.5 to 
2.8 A, inducing a transition from n-type Schottky con¬ 
tact to p-type Schottky contact at the interface. The 
conversion of Schottky contact type in G/P heterojunc¬ 
tion occurs when their interfacial distance increases to 
below 3.5 A. 

We find that when the interfacial distance increases 
larger than 4.5 A, the Dirac point of graphene is close 
to the conduction band of phosphorene, forming a n- 
type Schottky contact with a small barrier = 0.12 
eV at the interface as shown in Fig. 2d. That is be¬ 
cause the work function of graphene (4.3 eV) [59] is 
close to the nucleophilic potential of phosphorene (4.1 
eV). When the interfacial distance decreases from 4.5 
to 3.0 A, the charge transfer and chemical interactions 


between graphene and phosphorene are enhanced [113]. 
As the interfacial distance artificially decreases to 3.0 A, 
the Dirac point of graphene is close to the valence band 
of phosphorene, forming a p-type Schottky contact with 
a negligible p-type Schottky barrier of ^Bp = 0.06 eV at 
the interface as shown in Fig. 2d. Furthermore, strong 
interlayer-interactions even can induce weak n-type dop¬ 
ing in graphene when the interfacial distance artificially 
decreases to 2.8 A as shown in Fig. 2c. Therefore, G/P 
heteroj unction has tunable Schottky contacts and barri¬ 
ers as varying the the interfacial distance. 

Notice that the Schottky contact in G/P heterojunc¬ 
tion is very different to traditional metal-semiconductor 
Schottky ones in two important ways. One is that 
graphene is adsorbed physically on phosphorene in ul- 
trathin van der Waals heteroj unction without dangling 
bonds at the interface. Another is the Schottky con¬ 
tacts and barriers can be adjusted sensitively by varying 
the interfacial distance. Furthermore, atom doping in 
graphene and applying vertical electric fields can also be 
used to modify the work function of graphene and then 
to adjust the Schottky contacts and barriers in G/P het¬ 
erojunction. Therefore, G/P heterojunction can be used 
for tunable Schottky diodes in nanoelectronics. 
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g-CgN^ layer M0S2 sheet 


FIG. 5: (Color online) (a) Geometric structures and differential charge density of g-C 3 N 4 /MoS 2 heterojunction. The gray, 
blue, yellow and blue balls denote carbon, nitrogen, sulfur and molybdenum atoms, respectively. The red and blue regions 
indicate electron increase and decrease, respectively, (b) Density of states of g-C3N4/MoS2 heterojunction, g-C 3 N 4 and M 0 S 2 
monolayers. The Fermi level is set to zero and marked by green dotted lines, (c) Absorption coefficients of g-C 3 N 4 monolayer 
and g-C3N4/MoS2 heterojunction, (d) Schematic illustration of carrier transfer and separation in g-C3N4/MoS2 heterojunction. 
(Reproduced from Ref. [87], with permission from American Chemical Society.) 


Graphene/g-ZnO 

Monolayer g-ZnO is a wide-gap semiconductor. When 
forming a hybrid G/undoped g-ZnO heterojunction, the 
Fermi level still remains in the induced gap, showing that 
little charge transfer between graphene and undoped g- 
ZnO monolayer. This is also supported by the differential 
charge density as plotted in Fig. 3c. However, the in¬ 
homogeneous g-ZnO substrate can induce special charge 
redistribution in the graphene plane, forming intralayer 
electron-hole puddles as plotted in Fig. 3c. This effect 
can significantly enhance the electron conductivity and 
generate new photovoltaic activities [67]. 

For hybrid G/doped g-ZnO heterojunctions, we ob¬ 
serve that different charge rearrangement occurs in the 
heteroj unctions as shown in Fig. 3c. We find that A1 
and Li doped g-ZnO monolayers have a work function ei¬ 
ther 1.2 eV smaller or the same amount larger than that 


of graphene. Based on the Schottky-Mott model[lll], 
the Dirac points of hybrid G/doped g-ZnO heterojunc¬ 
tions are moved to below and above the Fermi level, 
respectively. As shown in Fig. 3c, electron-hole pairs 
are well separated at the interfaces of hybrid G/doped 
g-ZnO heteroj unctions. The calculated charge carrier 
density is about W(G/n-type g-ZnO) = 3.5x10^^ 
and Af/j(G/p-type g-ZnO) = 1.1x10^^ cm~‘^. Notice that 
these charge carrier concentrations are more than 3 or¬ 
ders of magnitude larger than the intrinsic charge car¬ 
rier concentration of graphene at room temperature (n 
= /6hiyp = 6x10^° cm~‘^). Therefore, the charge 

transfer in hybrid G/doped g-ZnO heteroj unctions can 
induce effective and tunable electron or hole doping in 
graphene for graphene-based optoelectronic devices, such 
as Schottky diodes and p-n junctions. 
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Graphene/MoS2 

High-performance field-effect tunneling transistors 
have been achieved experimentally [73, 74] in hybrid 
graphene and M 0 S 2 heterojunctions. Thus, the elec¬ 
tronic properties of hybrid G/M 0 S 2 and G/M 0 S 2 /G het¬ 
erojunctions affected by applying vertical electric fields 
are very desirable as shown in Fig. 4c. We find taht neg¬ 
ative vertical electric fields can induce p-type doping of 
graphene in hybrid G/M 0 S 2 heterojunctions. But, posi¬ 
tive electric fields almost have on effect on the electronic 
properties of hybrid G/M 0 S 2 heterojunctions. This is 
because electrons can easily from the Dirac point of 
graphene to the conduction band of M 0 S 2 but diffi¬ 
culty from the valence band of M 0 S 2 to the Dirac point 
of graphene due to the work function (4.3 eV) [68] of 
graphene close to the electronic affinity (4.2 eV) of mono- 
layer M 0 S 2 . 

We find that vertical electric fields can generate strong 
p-type but weak n-type doping of graphene at both 
negative and positive electric fields due to the symme¬ 
try in hybrid G/M 0 S 2 /G heterojunctions. The calcu¬ 
lated charge carrier concentrations in hybrid G/M 0 S 2 
and G/M 0 S 2 /G heterojunctions are shown in Fig. 4c. 
These values are more than 3 orders of magnitude 
larger than the intrinsic charge carrier concentration of 
graphene at room temperature (n = irkgT'^/dhvp = 
6x 10^*^ cm~‘^). The doping charge carrier concentrations 
of graphene in hybrid G/M 0 S 2 and G/M 0 S 2 /G hetero¬ 
junctions are increased with the vertical electric fields. 
Furthermore, the electron-hole pairs are well separated 
in hybrid G/M 0 S 2 /G sandwiched heteroj unctions with 
more excellent applications[73, 74] compared with hy¬ 
brid G/M 0 S 2 heterojunctions. Therefore, the field-effect 
in hybrid graphene and M 0 S 2 heterojunctions is effective 
and tunable for high-performance FETs and p-n junc¬ 
tions. 

g-C3N4/MoS2 

Monolayer g-G 3 N 4 and M 0 S 2 are all semiconductors 
with direct band gaps of 2.7 and 2.0 eV, respectively, g- 
G3N4 can be used to achieve both half reactions of water 
splitting under visible light, meaning that its bandgap 
cover both the water reduction and water oxidation po¬ 
tentials. However, the photocatalytic H 2 production ac¬ 
tivity of pristine g-G 3 N 4 remains poor, strongly depend¬ 
ing on the surface co-catalysts. On the other hand, M 0 S 2 
has been extensively investigated as a promising electro¬ 
catalyst for H 2 evolution due to its high abundance and 
low cost. According to some experimental studies, M 0 S 2 
can serve as a suitable co-catalysts to composite with 
g-G 3 N 4 to show enhanced photocatalytic performance. 

Our calculated band alignment between the g-C 3 N 4 
and M 0 S 2 monolayers reveals that the conduction band 


minimum and valence band maximum of the g-G 3 N 4 
monolayer are higher by about 0.83 eV and 0.15 eV re¬ 
spectively than those of the M 0 S 2 monolayer. This pre¬ 
dicted type-H band alignment ensures the photogener¬ 
ated electrons easily migrate from the g-C 3 N 4 monolayer 
to the M 0 S 2 monolayer, and leads to the high hydrogen- 
evolution reaction activity. The charge transfer between 
g-G 3 N 4 and M 0 S 2 results in a polarized field within the 
interface region, which will benefit the separation of pho¬ 
togenerated carriers. The calculated optical absorption 
curves verify that this proposed layered nanocomposite 
is a good light-harvesting semiconductor. These find¬ 
ings indicate that M 0 S 2 is a promising candidate as a 
non-noble metal co-catalyst for g-C 3 N 4 photocatalysts, 
and also provide useful information for understanding the 
observed enhanced photocatalytic mechanism in experi¬ 
ments. 


SUMMARY AND CONCLUSIONS 

In the present work, we review our recent theoret¬ 
ical studies on the structural, electronic, optical and 
electrical properties of 2D van der Waals heterojunc¬ 
tions via density functional theory calculations, in¬ 
cluding the Graphene/Silicene, Graphene/Phosphorene, 
Graphene/g-ZnO, Graphene/MoS 2 and g-G 3 N 4 /MoS 2 
heteroj unctions. We find that the weak interlayer van 
der Waals interactions in these 2D heteroj unctions can 
induce new properties and phenomena, such as bandgap 
opening, charge transfer and new optical absorption. 
Furthermore, our recently studied 2D van der Waals het¬ 
erojunctions have some special features which distinguish 
them from other previously studied 2D van der Waals 
heteroj unctions. With excellent structural, electronic, 
electrical and optical properties combined, 2D van der 
Waals heteroj unctions are expected to be with great ap¬ 
plications in efficient electronic, electrochemical, photo¬ 
voltaic, photoresponsive and memory devices. 
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APPENDIX A. SUPPLEMENTARY MATERIAL 

Supplementary material associated with this re¬ 
view article can be found, in the online version, at 
http://dx.doi.org/10.XXXX/XX.XX. 
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